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We study the possibility that dark energy decays in the future and the universe stops accelerating. 
The fact that the cosmological observations prefer an equation of state of dark energy smaller than 
-1 can be a signal that dark energy will decay in the future. This conclusion is based in interpreting 
a w < — 1 as a signal of dark energy interaction with another fluid. We determine the interaction 
through the cosmological data and extrapolate it into the future. The resulting energy density for 



dark energy becomes = a 



-3(X+«>d)_-/9(o-l) 



i.e. it has an exponential suppression for a ^> a — 1. 



In this scenario the universe ends up dominated by this other fluid, which could be matter, and the 
universe stops accelerating at some time in the near future. 



I. INTRODUCTION 

In the last few years the existence of dark energy as a 
fluid with negative pressure that accelerates the universe 
at present time has been established [1],[2]. Within the 
context of field theory and particle physics it is appeal- 
ing to interpret the dark energy as some kind of particles 
that interact with the particles of the standard model 
very weakly. The weakness of the interaction is required 
since dark energy particles have not been produced in the 
accelerator and because the dark energy has not decayed 
into lighter (e.g. massless) fields such as the photon. 
Perhaps the most appealing candidate for dark energy is 
that of a scalar field, quintessence [3], which can be ei- 
ther a fundamental particle or a composite particle [4] . It 
was common to assume the interaction between the dark 
energy and all other particles to be via gravity only, how- 
ever recently interacting dark energy models have been 
proposed [5]- [9]. The interesting effect of this interaction 
is two fold. On the one hand, the interaction between 
dark energy and matter, which can be for example dark 
matter or neutrinos [7] , is to give an apparent equation of 
state of dark energy smaller (more negative) than with- 
out the interaction and can be even smaller than -1 [9], 
as suggested by the cosmological observations. On the 
other hand it is also possible to have dark energy inter- 
acting with neutrinos and it is tempting to relate both 
energies since they are of the same order of magnitude 
[7] and a mass of neutrinos larger than 0.8eV imply that 
the dark energy cannot be a cosmological constant [8]. 

In general fluids with w < — 1 give many theoretically 
problems such as stability issues or wrong kinetic terms 
as phantom fields [10]. However, interacting dark energy 
is a very simple and attractive option which we will use 
in this letter. 

Since the dark energy dilutes slower than matter we ex- 
pect it to dominate the universe at late times. So, once 
the universe begins to accelerate due to dark energy we 
expect it to maintain this state of acceleration in the fu- 
ture and the universe will end up completely dominated 
by dark energy. In this letter we would like to study if 
this fate of the universe is unavoidable or we could have 



a transition from an accelerating universe to a non accel- 
erating one in which the dark energy decays into another 
fluid. We will show that the fact that the cosmological 
data [1], specially the SNla data [2], prefer an equation 
of state w of dark energy smaller than minus one can be 
a signal that dark energy will decay in the future and the 
universe will stop accelerating. This conclusion is based 
in interpreting ara<-lasa signal of dark energy inter- 
action with another fluid. We determine the interaction 
through the cosmological observations and extrapolating 
it into the future. 

This letter is organized as follows. In section II we 
present the generic evolution of two interacting fluids. 
In section III we introduce the effective and apparent 
equations of state and in section IV we give the evidence 
a dark energy decay. Finally we present our conclusions. 



II. INTERACTING DARK ENERGY 
A. Fluid Evolution 

The evolution of two interacting fluids p^ and pb, which 
can be quintessence scalar field <j> for dark energy "DE" 
and another fluid (pb), as for example matter or radia- 
tion, is given by 



p^ = -3Hp^(l + iity) - 5(t) 
p b = -3Hp b {l + w b ) + S(t) 



(1) 
(2) 



with H = a/a the Hubble parameter and S(t) the inter- 
action coupling. This S is a dissipative term, it depends 
on the interaction term between the particles of p^ and 
Pb and is in general a function of time. The equation 
of state parameters w$ = p^/p^ and Wb = Pb/pb may 
be functions of time. Without lack of generality we will 
take Wb > Wfj,, which is consistent with assuming p^ as 
the dark energy in the absence of an interaction term. 

A simple general solution for pb can be obtained by 
taking Wb constant and 5 = A(t)pb- Eq.(2) becomes pb — 
— [3H(1 + w#) — A]pb and it has a solution 



Pb 



p bo a-^ 1+ ^ Jl dtA 



Pb& 



dlA 



(3) 
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with pt = pb a~ 3 ( 1+Wb \ the evolution of pb without an 
interacting term, i.e. 5 = 0. We take a(t Q ) = a Q = 1 
at present time. A similar solution can be obtained for 
P0 with constant and setting 5 = Ap^ the solution to 
eq.(l) gives 



-m+^) e -IL dtA = n^e-C dtA 



(4) 



with p0 = p^a 3 ( 1+M, «) the evolution for 6 = 0. Of 

course A, A are related by 5 = Ap^ = Apt,. Clearly the 
sign of A, i.e. 5, determines whether p^ and pb will evolve 
faster or slower with the interaction term than without 
it. So, for A positive pb will dilute slower while for A 
negative it will dilute faster. 

If we take the ratio y = p^,/ pb = fl<p/flb from eqs.(3) 
and (4) we have 

y = P± = y a 3 ^-™*) e" SL *^(i+V») (5) 
Pb 

with y = p<f> /pbo, and taking the derivative of y w.r.t. 
time we find 



y = 3Hy [Aw - T] 

with Aw = Wb — Wrf, and 

T = J_ ( P4> + Pb \ = A(l + y) 
3H \ p^pb J 3Hy 



(6) 



(7) 



The value for y is constraint to < y < oo with y = for 
= and y = oo for pb = 0. Clearly from eq.(6) we see 
that the evolution of y depends on the sign of Aw — T. 



1. Non Interaction solution: Aw > T 

If (wb — w$) > T then y is positive and y will increase, 
i.e. pb will dilute faster than p^, and we will end up with 
Ptj, dominating the universe. In this case the interaction 
term 5 is subdominant and the evolution of p^ and pb is 
the usual one, i.e. p^ oc a^ 1+w ^ and pb oc a~ 3 ( 1+Wb K 



2. Interacting solution: Aw < T 

For pb to dominate the universe we need y -C 1 at 
late time and the (interaction) term T should dominate 
over Aw. A simple example is when A is constant and 
positive. In this case eq.(5) is 



a 

Vi 1 — 

a; 



e- M 



(8) 



3. Finite solution: Aw = T 

A solution to eq.(6) with y constant and pb ^ 0, p^ ^ 
is only possible if A/H (or A/H) is positive and constant 
since T = A(l + y)/3Hy = 1(1 + y)/3H = Aw must be 
constant and positive, taking Wb — w^ > constant. Let 
us take A/H = C > constant, i.e. with S = Ap^ = 
CHp^, and from y = 0, i.e. T = A(l + y)/SH = Aw, we 
get a stable value of y given by 



Vs = = T 3 ( Wb 

"6 A 



3(w b - w<f>) 
C 



-1. (9) 



It is easy to see that the solution y s is stable since from 
eq.(6) the fluctuation Sy = y — y s to first order behaves 
as dy/Sy = —Ay s (A is positive by hypothesis) giving 
5y — > 0. Furthermore, since A is proportional to H = a/a 
then we can integrate eq.(4) with Adt = [3(w(, — w ( f,)/(l + 
y s )]Hdt = CHdt = Cda/a to give 



P4> = P<i>% 



-3(l+«v)-C 



= P4». 



-3(1+ 



(10) 

where we have used that C = 3(wb — w^)/(l + y s ). Since 
C is positive the solution in eq.(10) dilutes faster than 
the non interacting solution ps oc a~ 3 ( 1+M, «) and with an 
effective equation of state w e ff = w^ + C/3 = (wb + 



III. EFFECTIVE AND APPARENT EQUATION 
OF STATE 

We will now introduce the effective cq. of state w e ff 
and the apparent eq. of state w app . 



A. Effective Equation of State 

To obtain an effective equation of state we simply 
rewrite eqs.(l) and (2) as 



/90 = SHp^l + Weff) 
pb = -3Hp b (l+Wbeff) 

with the effective equation of state defined by 
<5 S 



Weff =W 4> + ' 



Wbeff = w b 



3H Pb 



(11) 



(12) 



at late times for any value of Wb — w,p. 



3Hp A ' 
The solution to eq.(ll) is 

p = e -3 £(1+10.// )*./«. = a -3(l+^) e "£^ d * 

p b = e -3f* o ^+w beff )da/a = a -3 [l+Wb) Jlj- h dt ^ 

where we have assumed in the second equality of eqs.(13) 
that w,p 7 Wb are constant. We see from eqs.(13) that 
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w e ff,Wbeff give the complete evolution of and pi,. 
For 5 > we have w e ff > w^ and the fluid p^ will di- 
lute faster then without the interaction term (i.e. (5 = 0) 
while pb will dilute slower since Wbeff < Wb- Which fluid 
dominates at late time will depend on which effective 
equation of state is smaller. The difference in eqs.(12) is 

Aw e ff = Wbeff - W e ff = Aw - T (14) 

with T defined in eq.(7) while the sum gives 

QbWbeff + tt<t>Weff = tt b Wb + fi^Uty. (15) 

Clearly the relevant quantity to determine the relative 
growth is given by T and if T > Aw we have Aw e f / < 
and pb will dominate the universe at late times while for 
T < Aw we have Aw e ff > and p^ will prevail. For no 
interaction S = and T = giving Aw e ff = Aw > 
and p$ dominates at late times. If T = Aw then Wbeff — 
w e ff and the ratio of both fluids pb/p<j> will approach 
a constant value, and if the universe is dominated by 
+ pb, i.e. 0,$ + Qb = lj then eq.(15) gives 



< Weff = w b n b + w^(l - fi 6 ) < w b , 



(16) 



i.e. the effective equation of state is constraint between 
w<f, and Wb- Of course eqs.(14) and (7) are consistent with 
the analysis of eq.(6). 



B. Apparent Equation of State 

An interesting result of the interaction between dark 
energy with other particles is to change the apparent 
equation of state of dark energy [5]- [9]. An observer that 
supposes that DE has no interaction sees a different evo- 
lution of DE as an observer that takes into account for 
the interaction between DE and another fluid. This ef- 
fect allows to have an apparent equation of state w < — 1 
for the "non-interaction" DE [9] even though the true 
equation of state of DE is larger than -1. 

Let as take the energy density p ~ p^ + Pb = Pde + Pb- 
The energy densities p<f,,Pb are given by eqs.(l) and (2) 
and these two fluid interact via the 6 term. On the other 
hand the energy densities pde and p~b do not interact 
with each other by hypothesis and therefore we have pi = 
-3H(1 + w b ) and p DE = -3Hp DE (l + w ap ), i.e. 



Pb = Pb Q a 
Pde = PDE a 



3(l+u)(,) 
3(l+to ap ) 



(17) 



if Wb,w app are constant. It was pointed out that the 



apparent equation of state w app 
than -1 and it is given by [9] 



can take values smaller 



War, = 



x ^_p±(pb_ l 

Pb 



1-x 
Pb a 



(18) 



f* dtA -, 

jJto — 1 



sit) 


w e ff, Ma 


w apP (a < do) 


W a pp(a > do) 


5 > 
<5 < 


W ef f > Wj, 
W e ff < W4, 


W ap p < U>0 
Wapp > «I0 


W a pp > W^y 
Wapp < U>0 



TABLE I: We show the different relative sizes of w app and 
Wef f with respect to w<t, as a function of the sign of the inter- 
action term S(t). 



valid if Wb — and we have used eq.(3). We see from 
eq.(18) that for pb < p~b, i.e. A > (or 5 > 0), we have 
x > and w ap < w^ for t < t which allows to have a 
w ap smaller than -1. 

The result in eq.(18) can be generalized to the inter- 
action between two arbitrary fluids. Taking the time 
derivative of p = p^ + pb and using eqs.(l) and (2), for an 
arbitrary interacting term S, we get pde = —'AH pde(1 + 



Wr. 



,) with 



w 



WbX 



ap 



1-X 



(19) 



with 



x = 



Pb 



Pb 
Pb 



Pb a 



-3(l+w b ) 



r 



dtA 



- 1 



(20) 

where we have used eq.(3). Of course w ap , x in eq.(19) 
reduce for Wb = to w ap , x as given in eq.(18). For t = t a 
we have x = and if there is no interaction 5 = A = we 
have pb = Pb, x = and w ap = w^. An apparent equation 
of state w ap < W0 is given for < xAw — x(w b — w$), 
i.e. x > 0. A positive x needs pi, < Pb which from eq.(3) 
implies a positive interaction term 5 = Apb- We also see 
that a positive Wb gives a more negative w ap than for 
w b = 0. 

It is interesting to note, see table I that the apparent 
equation of state w app is smaller than w^, for a positive S 
while the effective equation of state w e ff is in this case 
larger than w^. This clearly shows that the apparent 
equation of state is an "optical" effect not a true evolu- 
tion. 



IV. EVIDENCE FOR DARK ENERGY DECAY 

The SNla observations prefer an equation of state 
w < — 1 for dark energy. In principal a w < — 1 for a 
fluid is troublesome since it has instabilities and causal- 
ity problems. However, as seen in section III this can 
be an optical effect due to the interaction between dark 
energy with other particles as for example dark matter 
or neutrinos. 

Here, we will assume that w < — 1 due to this inter- 
action and we will show that this can be interpreted as 
a signal for a dark energy decay in the future with the 
universe no longer accelerating. 

The complete evolution of dark energy is given by the 
effective equation of state w e ff given by eq.(12) while the 
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observed equation of state, if we assume no interaction, 
is given by eq.(18), i.e. w app = - x). 

For values of the scale factor a close to present day 
a = 1 we propose to approximate x linearly by 



x = j3{a a — a) = —(3Sa 



(21) 



with (3 a constant to be determined by observations and 
Sa = a — a . A positive x in the past (i.e. a < a ) requires 
P > 0. Eq.(21) satisfies the requirements x(a D ) = and 



we have for 



J app 



> x < for a > a Q . 



The SNla data are in the range 1 > a > 2/5, i.e. for 
a redshift < z < 1.5, and the best fit solution has an 
average equation of state < w —1.1 [1]. Taking the 
average of w app = w^/ {1 — x) we have 



ippda w^Log [1 — /3(1 — <zi)] 
da /3(ai — 1) 



(22) 



where we have used eq.(21) and a a = 1. As an example 



let as take < 



■app 



>= —1.1, as suggested by the obser- 



vations [2], [8], and = — 0.9,Wh = 0, ai = 2/5. In 
this case we obtained from eq.(22) the value (3 — 0.56. 
If instead of taking the average as in eq.(22) we simple 
take the w app evaluated at the extreme points we have 



^ Wapp ^ 2^^ a PP^°^ ^app 

(ai)) giving the simple 
analytic expression (3 = {1 _ ai )(2< Wapp >~ w ,) 



and for our 

previous example we find (3 = 0.51. The difference in 
determining (3 is small and the expression for [3 becomes 
analytic as a function of < w app >. 

Now, we would like to determine the effective (true) 
eq. of state for dark energy. From eq.(18) we have 



-A8t 



(23) 



where we have approximated e^*° dtA ~ 1 — A5t and St = 
t — t . Taking S = Apt, and eq.(23) we have S — —xp^/dt 
and using St = 5a/ aH and eq.(21) we get an interaction 
term 



6 = af3Hp c / } . 



Eq.(12) becomes then 



Weff =W4,+ 



5 a(3 



(24) 



(25) 



and the effective equation of state for the b-fluid gives 

5 a/3 

Wbeff = Wb -mp- b ^ Wb -YYh- (26) 

Using eq.(4) or equivalently eq.(13) with the interaction 
term given in eq.(24) we get an energy density 



H = a -3(l+^) e -/3 (a -l) 



(27) 



which shows that p^ dilutes as a 3 ( 1 + w <i>) for a <C a Q — 1 
and pd, is exponentially suppressed for a 3> 1. From 



eqs.(25) and (26) we find that Wbeff < w e ff for a > 
a* = 3AwO^//3 (taking fl^ + = 1), i.e. for a > a* 
dark energy dilutes faster than the b-fluid. 

For our previous example with = — 0.9, w;, = 0,/3 = 
0.56 and present day values f^ = 0.7 and f2{, = 0.3, we 
find w e ff = Wbeff at a* ~ 1.3. Furthermore u> e // = at 
a = — iw^l (3 = 4.8, i.e. due to the interaction w e // grows 
from w e ff — —0.9 to w e ff — 0. In this case the universe 
goes from a decelerating to an accelerating epoch and 
back to a decelerating one at a scale factor = 5.1. 

From eq.(26) with w b — we have Wbeff < 0. Since 
Wbef f < w e f f for 00 ^ and a > a* then we will have 
at late times — > and Wbeff — * with a universe 
completely dominated by the fluid pb, matter in this case, 
and decelerating. 

We show in fig.l the evolution of £1$, fib in the example 
with an interaction term 5 = —xp^/St = —af3Hp c p with 
x given by eq.(21) and = —0.9 and (3 — 0.56, such 
that the average < w app >= —1.1. Notice that for a < 1 
the dark energy grows relative to the b-fluid while for 
a > 1 (i.e. in the future) the b-fluid dominates. In fig. 2 
we show the effective equations of state given by eqs.(25) 
and (26) and both are larger than -1 at all times. In fig. 3 
we show the apparent eq. of state (c.f. eq.(18)) where 
w a pp is smaller than for a < 1 as suggested by the 
SNla data but it becomes larger than Wcf, for a > 1. 
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FIG. 1: We show the energy densities fl^jQb (red (dotted) 
and blue (solid) respectively) as a function of the scale factor 
a. The fe-fluid dominates at late time. The vertical line is 
present time a — 1. 



V. CONCLUSIONS 

The cosmological observations prefer an equation of 
state w < — 1 for dark energy. We obtain an appar- 
ent equation of state for dark energy smaller than mi- 
nus one due to the interaction between dark energy and 
another fluid (b-fluid). Form the observational data we 
determine the interaction term, close to present day, 
and we show that this interaction imply that dark en- 
ergy will dilute faster than the b-fluid. The interac- 
tion term is S = afiHp^, which gives an energy density 
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o 

S -0.2 
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-0.4 
i -0.6 
-0.8 



2 3 
a 



-3(1+-!^) ~/3(a-l) 



it has an exponential suppres- 



sion for a a . The resulting universe is a decelerating 
universe dominated by the b-fluid, which could be dark 
matter. 



FIG. 2: We show the effective eqs. of state w e ff,Wb e ff (red 
(dotted) and blue (solid) respectively) as a function of the 
scale factor a. The vertical line is present time a = 1. 
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FIG. 3: We show the apparent eq. of state w apP as a function This work was also supported in part by CONACYT 

of the scale factor a. It is less than -1 for a < a = 1. The pro ject 45178-F and DGAPA, UNAM project IN114903- 
vertical line is present time a — 1. o 
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